Local properties of extended self-similarity in 3D turbulence 
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Using a generalization of extended self-similarity we have studied local scaling properties of incom- 
pressible homogeneous isotropic 3D turbulence in a direct numerical simulation. We have found that 
these properties are consistent with lognormal-like behavior of the velocity increments with moder- 
ate amplitudes for space scales r beginning from Kolmogorov length n up to the largest scales, and 
in the whole range of the Reynolds numbers: 50 < R\ < 459. The locally determined intermittency 
exponent fi(r) varies with r; it has a maximum at scale r = 1477, independent of 7?^. 

PACS numbers: 47.27.Ak, 47.27.Jv, 47.27.Gs, 05.20.Jj 



So-called extended self-similarity (ESS) in incompress- 
ible turbulent flows is intensively studied in recent years 
(see, for instance, and references therein). The ESS 
implies scaling relation between moments of different or- 
der. For example, for absolute value of longitudinal ve- 
locity increments over separation r in the inertial range 
the ESS means 



.p\p(q) / p(p) 



(1) 



where the scaling exponent p(q) is some function of q. It 
is shown in numerous experiments and numerical simu- 
lations that the range of applicability of the ESS is sub- 
stantially larger than that for ordinary self-similarity, and 
the ESS can exist even for situations where the ordinary 
self-similarity cannot be observed at all. 

The aim of the present Letter is twofold. The first is 
to study local properties of ESS defined as 



«) ~ < 



3\Cp(r) 



(2) 



where Cp( r ) are depending on r. The local approach al- 
lows to develop some old ideas. At past a representation 
Cp = p/3 — l-i>p(p — 3)/18 derived from a lognormal model, 
was shown to hold for the values averaged over an iner- 
tial interval of certain extension. Here we will show that 
the representation holds locally (i.e. with £ p (r) and /_t(r) 
depending on 7') in the above defined sense, ranging from 
the dissipative Kolmogorov length 77 to integral scale for 
any observed Reynolds number in direct numerical sim- 
ulations up to R\ = 459. 

The second aim is to examine the implication of the 
variation of the local intermittency exponent /i(r). It is 
often mentioned that there is no characteristic length in 
turbulence, so that the structure functions obey a power 
law in the inertial region and the associated scaling expo- 
nents are independent of scale. However, a recent study 
of turbulence reveals the existence of the structures [^). 
A question, then, naturally arises as to whether there 
may be any characteristic length ascribed to the struc- 
tures [fni. The peculiar variation of /_t(r) with respect to 



r, observed in the Letter, indicates that there is a certain 
length affecting the ESS. 



The She-Leveque model |ll) is also very popular in the 
last years. For this model the local exponent is expressed 
in a general form: 



Cp(r) = | (1 - 7) + _^_ 



P/3 



(3) 



where f3 and 7 can be also supposed to depend on r; /3(r) 
and j(r) could be in principle evaluated by plotting the 
local exponent Cp( r ) against p with r fixed. Since here 
we focus the lower order structure functions which are 
computed more reliably than the higher order ones, the 
comparison of the data with the formula is done using 
a simpler lognormal model with one parameter. If one 
could compute the higher order structure functions, the 
comparison should be made with the She-Leveque model. 

Let us begin with the derivation of useful formulae on 
the basis of the lognormal distribution of the dissipation 
rate e r averaged over spheres of radius r: 



P(s r ) = 



V27RT 



: exp 



(lne r — a)' 
2^2 



from which we obtain 



(e r ) q 



= <t 2 9 (<j-1)/2 



(4) 



(5) 



that results in a parameter-independent type of ESS 
of turbulent energy dissipation 



g(g-i) 
tP\ \ p(p-i) 



(6) 



What is the equivalent relation for the velocity incre- 
ment? According to the refined similarity method JT3 nfl] 
u r ~ (re,,) 1 / 3 in the inertial region, while u r ~ re\ in 
the dissipative region. If we take regions of scale r be- 
tween the inertial and dissipative regions, some are de- 
scribed by the former relation, while the others are by 



1 



the latter one, so that u r is related to e r in a probabilis- 
tic way. For the sake of simplicity we assume a general 
relation 



c(r)e 



l/a(r) 



(7) 



in a mean sense. Here a(r) is a function of r and coeffi- 
cient c(r) does not necessarily scale with r. Substituting 
(0) into (|) yields 



pip— a) 



xq/a 
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a\P/ a 
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Note that this expresses how the structure functions of 
different order are related to each other with value of r 
fixed. This relation can be considered as a functional 
equation and a solution to this equation is 



q/a(r)+b(r)q(q—a{r)) 



(9) 



where b(r) is an arbitrary function of r. Relation (^|) is a 
generalization of ordinary ESS 0,0- The difference be- 
tween ordinary ESS and relation (0) is that parameters 
used in (^|) can depend on r. Therefore we will call this 
type of ESS as extended local self-similarity (ELSS). 

To be consistent with the present data processing we 
express the gth order structure function in terms of 3rd 
order structure function as (||). Making use of (||), we 
are led to 



Cq(r) = I 



18 



Q(Q - 3), 



where 



H(r) 



6b(r)a(r) 



l + a(r)(3 - a(r))b(r) 



(10) 



(11) 



In the frame of ELSS the exponent ( q depends also on r, 
and below we compare (||) and (|Io| ) with data of DNS for 
different values of r. It should be noted that the ELSS 
expression (10) holds for any value of a(r), so that the 



expression can be compared with the data for any scale 
separation without paying attention to which region is 
being considered. 

We have performed a series of direct numerical simula- 
tions (DNS's) of incompressible homogeneous isotropic 
turbulence using a resolution up to 1024 3 . Reynolds 
numbers range from 50 to 459 jlq] . The random force is 
statistically homogeneous, isotropic and Gaussian white, 
and applied to the band l<fc<3 in which the forcing spec- 
trum is constant. The code uses the pseudo spectral 
method and the 4th order Runge-Kutta-Gill one. Ini- 
tial conditions are Gaussian random velocity fields with 
the energy spectrum E{k) oc k 4 exp(— 2(k/ko) 2 ), and the 
resolution is N = 256 3 for R\ = 69, N — 512 3 for 
Rx = 125,176,259, and N = 1024 3 for R x = 374,459. 



After about two eddy turnover times all the turbulent 
fields attained statistically steady states, which were con- 
firmed by observing the time evolution of the total energy 
and enstrophy, and the skewness of the longitudinal ve- 
locity derivative. For R\ = 459 run, the Reynolds num- 
ber was gradually increased through two steady states. 
The condition k max rj > 1 for the resolution of DNS is 
satisfied for most runs, but that of R\ = 459 is slightly 
less than unity. The statistical averages were taken as 
the time average over tens of turnover times for lower 
Reynolds numbers and over a few turnover times for the 
higher Reynolds numbers; over 45 samples during 2.9 
eddy turnover times for R\ = 374 and over 31 samples 
during 1.4 eddy turnover times R\ — 459. Computations 
with R\ < 259 have been done using a vector parallel ma- 
chine with 16 processors, Fujitsu VPP700E at RIKEN, 
and those for higher R\, using Fujitsu VPP5000/56 with 
32 processors at Nagoya University Computation Center. 

Now turn to the data analysis. Figure 1 is a plot of 
(ii 2 )/(e?7) 2 / 3 against rjr\ for various values of Reynolds 
number, where e is the average dissipation rate, and 
rj = (v 3 /e) 1 / 4 with molecular viscosity v. Here a straight 
solid line proportional to r 2 / 3 is inserted. It is remark- 
able that all data points collapse on a single line in the 
dissipative region, which indicates that all simulations 
are carried out with the good resolution at small scales. 
Although the slope of (u 2 ) could be estimated for large 
Reynolds numbers as seen from Fig.l, the scaling expo- 
nents of higher order structure functions as well as low 
order ones for small Reynolds numbers can be reliably 
evaluated only on the basis of the ESS method, i.e. by 
plotting ( U p) against (u 3 ) 

In order to know the r-dependence of £ p (r) for various 
Reynolds numbers, we prepare Fig. 2, in which C P (r) with 
p = 4, 6, 8 are depicted for R x = 69, 125, 259, 374 and 459. 
(The eighth order structure function is confirmed to con- 
verge statistically.) Note that the data for scales larger 
than integral scales are not shown, because a universal 
property of turbulence is not expected in those data. It 
is remarkable that there is a dip at about rjr\ ~ 10, and 
that it grows in depth with Reynolds number. The excep- 
tion is the case R\ = 69, where a dip does not appear. As 
the scale increases beyond the dip, C4( r ) an( i CeCO tend 
to approach constant values, although the corresponding 
data for R\ = 259 behave in a slightly different way from 
other cases. For p = 8 the situation is the same as for 
p — 6, but the variation is larger. It is of interest to notice 
that for the largest Reynolds number 459 the flat region 
is observed in the interval 100<r/?7<300, which may be 
identified with the inertial region. For smaller Reynolds 
numbers it is a little difficult to find the flat region. For 
R\ = 69 we see the flat region in the interval 10<r/?7<30, 
and the corresponding slopes for p > 4 are larger than 
those for R\ — 459. However, the flat region at smaller 
scales for R\ = 69 is different from one at larger scales 
for R x = 459. 
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Before we compare the formula (10) with the data, 
the validity of the lognormal distribution of the veloc- 
ity increment should be ensured. The pdf of the in- 
termediate amplitudes is certainly lognormal. In order 
to determine the range of the lognormality, we calcu- 
lated a peak position u*(p) of uP.P(u r ) at rjr\ = 24 for 
R\ = 121. The pdf is satisfactorily fitted by a log- 
normal curve for u*{p = 1.5) < u r < u*(p — 6). On 
the other hand, the pdf of e r is found to be lognormal 
in much wider interval. If we employ the same nota- 
tion as above, the lognormality of e r holds at least in 
the interval of e*{p = -4) < e r < e*(p — 6) for the 
same Reynolds number. If the refined similarity hy- 
pothesis vlf. ~ re r is assumed to hold for any ampli- 
tude, the corresponding pdf of u r should be lognormal 
for u*(p = —12) < u r < u*(p = 18), which is much wider 
than the observed lognormal interval. The reason for the 
discrepancy is that the refined similarity holds only for 
the intermediate amplitudes of e r in agreement with the 
observation jl5] |l7). Hence, the use of the lognormal ex- 
pression for the exponent (10) is completely justified for 
intermediate values of p. 

In order to analyze a nature of the r-dependence of the 
local slope Cp( r ), we rewrite ( |To|) in the following form: 

Hs^)-^ (12) 

What is the range of pi To decide the range we cal- 
culated C P (r) for various values of p at units of 0.1 at 
r/r) = 24 for R\ = 121, and plotted Cp/p against p. Al- 
though such a plot is not given here, Cp/p is fitted by a 
straight line for 1 < p < 6. The curve is deviated upward 
from the straight fitting line for p > 7 and downward for 
< p < 1. Hence the comparison will be made in the 
range 1 < p < 6, but the intermittency coefficient /i(r) 
estimated there plays a significant role to represent even 
the whole intermittency effect. 

Figure 3a shows a curve C, P {r)/p vs. p for several val- 
ues of r/r/ with R\ = 69. Straight lines in this figure 
are the best fit lines for 1 < p < 6. The data points for 
p = 7, which are not included for the comparison, are 
slightly deviated from the lognormal lines as mentioned 
above. Intermittency index, /i(r), can be calculated from 
this figure using the slope of the fitting straight lines or 
the intersection point of the fitting straight lines with 
vertical axis (cf.([l2])). The calculated values of n(r) are 
shown in the inset of the Fig. 3a (/ii corresponds to cal- 
culations using the 'slope' method, whereas \ii are given 
by the 'intersection' method.) In the plot of /i(r) Tay- 
lor microscale A and an integral scale L are marked for 
convenience. (We have confirmed the well-known pre- 
diction H A/ry = l^l^R 1 ^ 2 and L/rj ~ R 3 / 2 .) For 
Reynolds numbers 259 and 459 we have obtained similar 
pictures ( Figs. 3b and 3c). Other Reynolds numbers give 
the same result. 



As seen from Figs. 3a to 3c, fi(r) substantially de- 
pends on r with a typical "two-maxima" shape. Even for 
R\ = 69 in Fig. 3a there is a small peak around r/r] ~ 10. 
Local maximum of the \x(r) at smaller scales exhibits a 
few interesting properties. (I) For R\ — 50 the maximum 
is not observed, and the first appearance of the maximum 
occurs at R\ between 50 and 69. (II) A position of this 
maximum (normalized by rf) is independent of Reynolds 
number, and takes value £\ ~ 14?y, which is actually the 
beginning of the dissipative region. (Ill) Value of /i(^i) 
scales with R x as fi(h) ~ R° x 26 - ( iv ) The scale A is 
located on the right descending hill. 

The flat region where /j,(r) is constant appears in be- 
tween r/r] ~ 100 and rjr\ ~ 300 for the highest Reynolds 
number R\ = 459 (see Fig. 3c). In such an interval 
\i(r) w 0.25. This value is consistent with those known 
in literature for observations corresponding to very large 
Reynolds number pp| , pl| . 

It should be noted that strong tube-like vortices are 
believed to be sparsely distributed in space in fully de- 
veloped turbulence |^,^2|. Those vortices are frequently 
assumed as Burgers' vortex with mean radius IO77 [ p2[ . 
The energy dissipation takes place strongly around those 
vortices. Therefore extreme intermittency of the energy 
dissipation at this scale is consistent with the observation 
of the maximum of /i(r) at scale r = \^r\. On the other 
hand, the usual scaling region, where Cp( r ) as well as /i(r) 
are independent of scale in a certain interval of scale, can 
be seen only for largest Reynolds number, R\ = 459 in 
the present work. Such a scaling region is located for 
r > A; in our simulation the scaling region starts from 
approximately 2.5A. 
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FIG. 1. A plot of (v,r) divided by (en) 2 ^ 3 vs. r/n for var- 
ious values of Reynolds number. An inserted solid line is 
proportional to r 2 ^ 3 . Notice that all data points collapse on 
a single curve in the dissipative region. 
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FIG. 3. The ELSS exponents Cp(. r )/P against p obtained 
in the DNS for different values of r and for different Reynolds 
numbers; (a) R x = 69, (b) R x = 259 and (c) R x = 459. 
Straight lines, which are the best fit line for 1 < p < 6, in- 
dicate agreement of the data with the representation (12). 
The inset shows local intermittency exponent /i(r) calculated 
using the data. /Ui(r) corresponds to calculations using the 
'slope' method (described in the text), and fJ,i{r) does to those 
using the 'intersection' method. Calculated Taylor length A 
and the integral length L are marked for convenience. 
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